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Abstract 



The main problem in theoretical analysis of structures with strong 
confinement is the fact that standard mathematical tools: differential 
equations and Fourier's transformations are no longer applicable. In this 
paper we have demonstrated that method of Green's functions can be 
successfully used on low-dimension crystal samples, as a consequence of 
quantum size effects. We can illustrate modified model through the prime 
cubic structure molecular crystal: bulk and ultrathin film. Our analysis 
starts with standard exciton Hamiltonian with definition of commutative 
Green's function and equation of motion. We have presented detailed 
procedure of calculations of Green's functions, and further dispersion law, 
distribution of states and relative permittivity for bulk samples. After 
this, we have followed the same procedures for obtaining the properties of 
excitons in ultra-thin films. The results have been presented graphically. 
Besides modified method of Green's functions we have shown that the 
exciton energy spectrum is discrete in film structures (with number of 
energy levels equal to the number of atomic planes of the film) . Compared 
to the bulk structures, with continual absorption zone, in film structures 
exist resonant absorption peaks. With increased film thickness differences 
between bulk and film vanish. 



1 Introduction 



Interest in exciton subsystem studies appeared due to the fact that excitons 
are responsible for dielectric, optical (absorption, dispersion of light, lumi- 
nescence), photoelectric and other properties of crystals [1—3]. Studies of 
excitons in crystalline subsystems culminated with laser invention. 

In recent years theoretical investigations of quasi-two-dimensional exci- 
ton subsystems (nanostructures) were intensified, especially in the field of 
thin films, not only to obtain fundamental information regarding dielectric 
properties of these materials but also because of their wide practical use (na- 
noelectronics, optoelectronic [4—6], light energy conversion [9,10]...). What 
is unique for these structures is that they have changed properties compared 
to their bulk analogues [6—8]. 

We studied the basic physical characteristics of ultrathin dielectrics — 
molecular crystalline nanofilms [11,12], which could be used as surface layers 
for protection of electronic components or as special light filters. 
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This paper analyzes the influence of border-film structure on the energy 
spectrum of excitons (exciton dispersion law). Special attention was paid 
to the presence and spatial distribution of localized exciton states. Optical 
properties of these dielectric films were also investigated (their dielectric 
permittivity was determined). Results obtained in this work were compared 
with the similar results for the case of ideal infinite crystals, in order to find 
most important differences between these two systems. 

These analyzes may be conducted using methods of two-time, temperature- 
dependant Green's functions that are often used in quantum theory of solid 
state [13—15]. With adequately incorporated statistics, this method is being 
successfully applied in calculations of microscopic and macroscopic, as well 
as balanced and non-balanced properties of crystals. 

A question which justifiably arises is related to the mean of calculating 
Green's functions, which are "borrowed" from the quantum field theory and 
whose definition, i.e. usage, is based on variables with continuous spectra 
in unlimited (both direct and impulse) space! 

This work proves that the method of Green's functions may be successfully 
applied onto crystalline samples of such small dimensions that the quantum 
size effects are relevant [16]. In order to illustrate adaptation of this method, 
we will observe molecule crystal with a simple cubic structure: spatially 
unlimited (bulk) and strongly limited along one axis (ultra thin film). Our 
intention is to show a technique of application of Green's function onto 
spatially limited systems, so we excluded from calculation all really existing 
parameters: more complex crystalline structure, changes in boundary film 
parameters etc. 

2 Excitons in bulk-structures 

The discussion of dielectric properties of an ideal (with no defects, vacancies, 
etc) unlimited molecular crystal will start with standard exciton Hamilton- 
ian, which has following form in configuration space [1,15,17]: 

(1) H = H + Y J An P^Pn + X ™ P t P ™ + Z) Y ™ P t P n P t P rh , 

ft n,m n,fn 

where Pt and Pft represent creation and annihilation exciton operators at 
the node (site) n of the crystalline lattice. Quantity A^j represents the energy 
of the exciton localized at the n node, while the quantities X^^i and Yfi,m 
represent matrix elements of the exciton transfer from the node n to the 
node rh. 
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Properties of the model exciton system may be analyzed using commuta- 
tion Paulian Green's function [13,14,18]: 

(2) r m (t) = ((p H (t) i p+(o))) = e(t) <[p a (t), Pi(o)]) , 

which satisfies following equation of motion: 

(3) m|r^(t) = in<5(t)([p f j(t), P±(0)])+ 

+ e(i)([p 3 (i), ^] p+(o) - p+(o) [p*(t), ^]) . 

Using commutation relations for Pauli operators [18,19]: 

(4) [Pft,Pi} = {l-2PfPn)S nih ; 
[P H ,P A ] = [P n ,P+)=0; P| = (P+) 2 = 0, 

we obtain the equation of motion for Paulian Green's function: 

ihjT nih {t) = ihS(t)6 n rn (1 - 2<i£Pa» + 

(5) +A n T itrn (t) + J2Xra T TJ t )- 

i 

i i 

expressed by T mfih (t) = ({P+ {t)P R {t)Pft) | Pt(0)))- Paulian Green's func- 
tion of the higher (third) order. 

The basic problem with exciton theory is the fact that Pauli-operators P + 
and P are not Bose or Fermi operators, but a certain hybrid of both with 
a kinematics described by expression (4), that is Fermian for one mode and 
Bosonian for different modes. For precise analysis of exciton systems, which 
encompass effects of inter-exciton interaction, simple replacing of Pauli- 
operators with Bose-operators is not enough. Therefore, in Hamiltonian 
(1), Pauli-operators are replaced by their exact Bosonian represents [20]: 



(1 + z/)! 



2 

B 



(6) P + = B 



^ i 



(-2) 1 



oc 



P = y± 



[B Y B l 



' S 2) \. (B + Y +1 B v+l 



Our goal is to adapt Green's function method to spatially quantum (dis- 
crete, not continuous) structures and to see the influence of spatial limits 
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and disturbances of inner translational symmetry on changes of their macro- 
scopic physical properties. Paulian Green's functions from equation (7) will 
be therefore expressed using appropriate Bosonian Green's functions on the 
basis of approximate expressions following from (6): 

P^B-B + BB; P + ra B + - B + B + B ; 

(7) 

v ' P + P « B + B — B + B + BB . 

By this we obtain: 

r*a(t) = ((Pn(t) I Pt(0))) = ((B H (t) I B+(0)»- 

(8) - {(B H (t) | B+(0)B+(0)B^(0))) - «S+(t)Srt(t)SB(t) | B+(0)))+ 
+ «J3^(*)J3»(*)J3»(*) I B+(0)B+(0)B^(0))) . 

Further, by decoupling higher Green's functions using known Bose-commutation 
relations: 

(9) [B* Bf) = 6 ff ; [B U , B r ] = [b+, B+] = 
and by introducing retarded (Bosonian) Green's function: 

(10) ((B n (t)\B+(0)))=G m (t), 

terms in expression (8) become: 

((B n (t) | B+(0)B+(0)B*(0))) = Q^B^B+B+B^]) = 
= e(t) ({(S m + B+B H ) B+B M ) - (B+B+B rn B H )) = 2G Hrfl {t)M ; 

(11) ((B+(t)B n (t)B n (t) | B+(0)» = 2G r , A (t)No ; 
((B+(t)B n (t)B n (t) | B+(0)B+(0)B A (p))) = 2R HrA {t) G^(i) , 

where No is concentration of excitons, and Rftfft{t) is advanced Green's func- 
tion: 



No = {B + B) = ^J2{ e ^ )/e ~ 1 



(12) 

Rm,{t) = ((B+(t) | B A (0))). 

When expressions (10) and (11) are substitued into eq.(8) we obtain final 
expression for Paulian Green's function expressed using Bosonian Green's 
functions: 

(13) T HA (t) = (1 - 4M,) G m (t) + 2R HA {t) Gl A (t) + 0(N 2 ) . 

For Paulian Green's functions of higher order (7^-?-) at the left side of 
Green's function we simply replace Pauli operators by Bose-operators, and 
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on the right side approximation (7) takes place. In this way it follows: 

Taaic= i ^(o)» = m{t)Bz{t)B t (t) \ B$m- 

(14) - ((B+(t)B a (t)B^t) | B+(0)B+(0)B 5 (0))) = 

= N G G is (t) +M^G ss (t) - 2R S5 (t) G is {t) G dS {t) + 0«) . 

Expressions for T^, T-+~r~ and 7sv=t, which are expressed using Bosonian 
Green's functions, are substituted in a equation of movement for Paulian 
Green's function (5): 

\hj t [(1 - 4A/- ) G HrA (t) + 2R^{t) Gl~(t)] = md(t)6 n ri (l - 2{P+P H )) + 

[(1 - AM ) G m (t) + 2R Hrh {t) & m {tj\ + 

(15) +£^-[(1-^0) G r Jt) + 2R r ,(t)Gfjt) 

r 

-2^Xtf[Af Gfa(t)+M m G m (t) - 2R HA (t)G r Jt)G H , A (t)] + 
i 

+^ Y niWoG^(t)+M H ^ rn ,(t) - 2R t Jt)G Hni (t)G rn ,(t)] . 
i 

Since concentration of Frenkel's excitons in molecular crystals is very low 
(AT < 1 %), the equation above may be solved in the lowest approximation 
(this approximation is appropriate for neglecting anharmonic and nonlinear 
effects, i.e. non-calculating higher orders terms of exciton-exciton as well as 
exciton-phonon interactions, in agreement to estimates in [1—3,13—15]): 

(P+P n )n(B+B n )=M K0; 
A/^«0; G-G^O; G-R&O. 

"Decoupled" eq.(15) then obtains the following form: 

(16) \hjG HrA {t) = iM{t)5 nih + A H Gf; M (t) + Xnfilrn^) " 

I 



It is important to notice that this equation has to be obtained in the 
same form starting from effective Bosonian exciton Hamiltonian in harmonic 
approximation : 

n n,m 

and estimating Bosonian Green's function (10): 

GW(t) = ((Bft(t) | B+(0))) = 9(t) ([B n (t), B+(0)]), 
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with its equation motion: 

ih^G m (t) = \h8{t) ([B H (t), B+(0)])+ 



+ @(t)([B n (t),H ex ] B+(0) - B+(0) [B n (t),H ex ]) 



This (by time) differential equation (16) is solved using temporal Fourier's 
transformation: 

+00 +00 
(17) 4r(*) = / duf^u) e"^ ; 8(t) = ± [ e" 1 



2tt 



and thus we obtain 

oj) = 

2tt 



(18) hjjG nA {uj) = ^-Sftrfi + AftGfift(u)) + ^ X fU G lrn^) 



By using nearest neighbors approximation (I — >• it ± Aj): n ± Ai = n x ± 
1, n^, n z ; n ± A2 = n x , n y ± 1, n 2 ; n ± A3 = ra x , n^, n 2 ± 1 and taking into 
account that we are observing an ideal cubic structure, where exciton energy 
is the same at every node, and the energy transfer between neighbors is also 
the same: A^ = A; X ss± ^, = Xi i € {x,y,z}, equation above is taking 
following form: 

Gn x n y n z ,in x m y in z {^) = G rixriyriz ^jf i (uj) = 

= ^n I% « z ,™ + AG ni%nz ,,( W ) + 

(19) + X x [G nx+ i^ nynz .^(u) + G nx _i^ nynz .^{uj)\ + 

+ X y [G nxny+ i^ nz ^(uj) + Gn^-i^^ta)] + 
+ X z [G nxnynz+ i-^(uj) + G nxnynx -i^{u})\ . 

Since the crystal is unlimited, when solving this linear differential equation 
we may use the full spatial Fourier's transformation: 



By using these transformations and by composing the equation above, we 
obtain: 



(21) + 2 (Xj. cos a^fca; + Xy cos a y fcy + X z cos a z k z ) G^(co) 



huG^u) = - + AG £ (w)+ 
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and from there we may express Green's function: 
G-( w ) = — [Tiii) — A- 

(22) — 2 (X x cos a x A; x + X y cos a^A;,, + X 2 cos a 2 /c z )] _1 = 

in 1 
= 2^ frw-Ej; ' 

We obtain the energy spectrum in a bulk monomolecular crystal [13,14] by 
calculating real part of the pole of this Green's function: 

(23) Et = A + 2 (X x cos a^A^ + X y cos a^A^ + X z cos a z /c 2 ) . 

In order to perform comparison with dispersion law of excitons in the film, 
this expression will be written in a simpler (X x = X y = X z = — \X\, a x = 
o-y = a z = a) and non-dimensional form: 

E v - A 

(24) S % = = T xy + g z , 



where 



Fxy = — 2 (cos ak x + cos ak y ) ; Q z = — 2 cos ak z . 



This dispersion law is shown in Figure 1, as a function of two-dimensional 
value T xy : 

= &z {J~ xy ) , 

It is clear that for aki G [0, +ir], i = x,y,z (the first Brilouin zone), these 
values are within the intervals: 

^ G !" 4,+4 ! ; Uf S €[-6 I+ 6]. 
g 2 G [-2, +2], J 

The presence of permitted (continuous) energy levels is visible. 



Figure 1. Exciton dispersion law of bulk crystal 

Since molecular crystals are dielectric, it is essential to determine relative 
permittivity of these structures. The dynamic permittivity is defined by 
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systems response to external perturbation [1,15,20], by expression: 

ttS\X\ 



(25) 



-V) = i 



m 



;r( w ) + r(- w )] , 



where 5 is a frequency property of a given crystal and external variable 
electromagnetic field. As it was mentioned before, in this (zero) approxi- 
mation Paulian Green's functions (r) are transforming into Bosonian ones 
(G), therefore: 

(26) e- x (u) = 1 - ^ffl [G(u) + G(-lo)] . 

Substituting (22) with (23) in (26) and by denoting = |/| and — — = 

I I I I 

we obtain the expression for dynamic permittivity in a bulk sample of 

monomolecular crystal (molecular crystal with simple cell): 



(27) eAu) 



l + S\X\ 



E T 



(M 5 



El 

k 



-1 



1 + 5 



|/| 2 -(%+N)' 



Dependence of relative dynamic permittivity, expression (27), on reduced 
frequency (non-dimensional factor Hlj/A) of the external electromagnetic 
field is shown in Figure 2. 

The presence of a single absorption zone is visible within certain boundary 
frequencies. This energy zone was calculated for two-dimension center of 
Brilouin's zone (k x = k y = 0; k z = [0, ir]). For all other energies this crystal 
is transparent and has no spatial non-homogeneousness. 




r44 46 



Figure 2. Relative permittivity of bulk crystal 
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3 Excitons in thin film-structures 



Opposite to ideal unlimited structures, real crystals have no translation 
invariance property. The presence of certain boundary conditions is one of 
reasons for symmetry breaking [4—7]. Lets observe the ideal ultrathin film 
with simple cubic structure, made in the substrate i.e. by doping process. 
Here, by the term "ideal", we wanted to denote that there exist no breaking 
of inner crystal structure (no defects, ingredients, etc), and not as in no 
spatial limits. The film dimensions are such that it is unlimited in XY 
planes, and in z-axis has final thickness (L). This means that this film 
has two unlimited boundary planes parallel to XY planes, for z = and 
z = L = n z . 



3.1 The model 

The film-structure with primitive crystalline lattice (one molecule per ele- 
mentary cell): monomolecular crystalline film, with specified parameters is 
shown in Figure 3. 

Since boundary planes of the film are taken as being normal to z-axis, 
the index of parallel XY planes n z has values n z = 0,1,2, ... , N where 
N £ [2, 8] is for ultra thin films. Indices n x and n y , determining position 
of a molecule in every XY plane may have arbitrary whole-number values 
(practically from — oo to +oo). 



■ 




A 


n,=N 






X z 


A 


n z =N-\ 






X 


A_ _ 




L 




x z 










x z 


"a - " 


ii- =1 






x z 


A 


ii : =0 







A 


X/Y 



Figure 3. Cross-section of crystalline film- model 

For calculating exciton energies in this film, we start from equation (19) 
on which, due to spatial limitations of the film at z-direction, we may apply 
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partial spatial Fourier transformation: 

r _ _ _ \ J.k x a x (n x -m x ) ik y a y (n v -my) r 
°nm — AT »t / , c e °n z m z 

(28) 

x y h h 

(only along x and y directions). For expression shortening , it is convenient 
to introduce markings G nzmz (k x ,k yi uj) = G nzrriz . For X x = X y = X z = 

- \X\, a x = a y = a z = a and JC = lvl ^n z ,m z , we obtain: 

Z7T \X | 

(29) G nz _i )?Ttz + qG ' n z m z G riz j r \^ rriz — KL8 n , z >mz , 

where denotation is introduced: 



+ 2 (cos ak x + cos a/c^ 



_ E k~ A 

(30) = r— J- xy = £j: — T xy 

(the quantity J- xy is defined in expression (24), while £? i.e. E? which express 
a possible energy of excitons in films, will be calculated latter). 

Equation (29) is in fact a system of N + 1 nonhomogeneous algebraic 
differential equations with (starting boundary) conditions: G nZ)TrLz = 0, for 
n z < and n z > N + 1: 

2 G + Gi = /C 
Go + q Gi + G2 = /Ci 
Gi + ^G 2 + G 3 = /C 2 



G Uz -2 + QG nz -i + G„ 3 = /Cn z -1 
(31) G nz -\ + ^G„ 3 + G nz+ i = K riz 

G Uz + QGn z +l + G nz+ 2 = K-n z +l 

Gn-3 + QGn-2 + Gn-1 = ICn-2 

Gn-2 + qGn-i + Gn = /Cjv-i 
Gat-i + qGn = K-n 

where is describe: G nz = G nzrriz and IC5 nzt m z = JC nz , because m z is blinding 
index. 
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3.2 Dispersion law 



In order to find exciton energies, we need poles of Green's functions, which 
are obtained when a determinant of a system (31) is equalized with zero, 
i.e. 

• 

• 



(32) 



Dn+i(q) 



2 1 

1 Q 1 
1 Q 1 













1 Q 1 

1 e 1 

1 Q 



N+l 



and this determinant is actually a note for Chebyshev's polynomials of sec- 
ond kind (and (N + l)-th order) [21,22]: 

(33) D N+1 (g) = C N+ i(g) = qC n (q) - C N -i{q) . 

The condition Dn+i{q) = is reduced to Cn+i{q) = and is satisfied by 
N+l solutions in form: 

V 7T 



(34) 



-2 cos ■ 



iV + 2 ' 

Using this and replacing equation (30) we find: 
(35) 



l,2,...,iV + l. 



Ek x k y ty) — A — 2\X\ I cos ak x + cos ak y + cos 



N + 2 



In order to compare with dispersion law of excitons in a bulk we will write 
this expression in more simple, non-dimensional form {g u = —2cosa z k z {[i)): 

(36) £k x k y {v) = F xy + Qziy) ; Q z {v) = - [1 - cos a z k z (v)\ . 



Previous expression represents the dispersion law of excitons of ideal 
monomolecular film and has the same form as the expression (24) obtained 
for corresponding ideal unlimited structures, with difference that in (24) k z 
is practically continuous variable (interval [0, 7r/a]) as k x and k y , while here 
is discrete and is given by expression: 



(37) 



k z (y) = - 



a N + 2 ' 



1,2, ... ,N + 1. 



Graphical representation of dispersion law is given on Figure 4, showing 
possible exciton energies in ideal five-layer monomolecular film (full lines) 
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together with bulk boundaries (dotted lines). Same as for corresponding 
bulk crystalline structures, we will use ordinate for values of reduced non- 
dimensional energies £ u , depending on two-dimensional function JF X y, for 
which is used graph abscissa. 




Figure 4. Exciton energies of molecular film 



These analyzes have shown significant differences regarding dispersion 
law for excitons in spatially strongly limited systems (nanofilm-structures) 
as strictly result of border presence of these structures, in which energy 
spectra is highly discrete and have two gaps. Sizes of gaps depends on film 
thickness and decrease rapidly with its increase. 



3.3 Spectral weights 



In order to find certain Green's function we will start from the system of 
equations (29), which is now suitable to represent in the operator form: 



(38) 



V 



n+ivn+i 



JCn 



+i 



Gn+i = ^V+i^Qv+i , 



where: £>at+i is a matrix corresponding to the determinant of system Djv+i, 
Gn+i and /Qv+i are vectors of Green's functions and Kronecker's delta sym- 
bols: 

Gl,m z 



(39) 



/N+l 



G 



n z ,m z 



JCn 



ih 



+1 



2n\X\ 



Si 



J n z ,m z 



\ GN,m z J \ $N,m z j 

Since the inverse matrix 2^+1 ma Y be expressed through adjunct matrix, 
whose elements Dn~ are cofactors of elements from the direct matrix, we 
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may write: 



(40) = - vl +1 2*x\ ^ Dn *« 5 * m - = 

ih Dn z ,m z 

2tt \X\ fjv+i 

Cofactor D nzTnz calculation is based on knowing the system determinant 
Dn+i- 

Since for the equilibrium processes within system are important only di- 
agonal Green's functions G nz - n , z = G Uz , calculating cofactors D nz ,m z = D n , 
is significantly simplified. It turns out that they are equal to the product of 
two auxiliary determinants: 



(41) D nz = C nz C 



N-n- 1 



where C nz -\ and Cn-u z ar e corresponding Chebyshev's polynomials of sec- 
ond kind. Therefore, Green's function of the ideal film is: 



(42) G nz ik C - CjV -- 



2vr|X| C 



N+l 



These Green's functions are multipolar, since denominator consists of a 
polynomial Cn+i of N + 1-th order. Therefore factorization on simple poles 
must be performed [21]: 

N+1 ( \ 

f A o\ n 9n z ;nASu) 



2'k\X\ ^— ' Q — Q v 



Spectral weights g nz ;n z (Qu) = 9n- then may be expressed using: 
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Using the rule for derivative of determinant we obtain: 



(45) 



T q Cn+i{q) 



N+l 

= ^2 Ci-l(Qv)CN+l-i(Qv 
Q=Qv i=l 



N+l 

£ 

i=i 



s'mi^ u sin(iV + 2 — 
sin 2 i v 



= sm 2 (n z + 1)6, 

N+l 

^2 sin2 

i=l 

and spectral weights became: 



N + 2 ' 



(46) 



JV + 2 



■ sin 



(n z + 1) 



N + 2 



The spectral weights of Green's functions are squares of the modules of the 
wave function of excitons [1—3,13—15] and enable determination of spatial 
distribution, i.e. probability to find excitons with certain energies per layers 
of crystalline film. This is in fact the spatial distribution of probability to 
find certain energy state of excitons. 

Numerically calculated, values of reduced energies and corresponding 
spectrum functions (spatial distribution of probability) for four-layered film 
(N = 4, where k x = k y = 0) are shown in the table. 

Table 1. shows spatial distribution of exciton energies occurrence proba- 
bilities in ideal monomolecular film. 



Reduced 
relative 

ENERGY 


ULTRATHIN FILM 
atomic plane 


upper 
boundary 


first 
inner 


central 
(inner) 


.last 
inner 


lower 
boundary 


-1,73205 


0,08333 


0,25000 


0,33333 


0,25000 


0,08333 


-1,00000 


0,25000 


0,25000 


0,00000 


0,25000 


0,25000 


0,00000 


0,33333 


0,00000 


0,33333 


0,00000 


0,33333 


1,00000 


0,25000 


0,25000 


0,00000 


0,25000 


0,25000 


1,73205 


0,08333 


0,25000 


0,33333 


0,25000 


0,08333 



Table 1. Exciton probabilities in the ideal four-layered film 



This table shows that for one certain energy, probability of exciton oc- 
currence per all layers is equal to one, and that probability per one layer for 
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all energies is also equal to one, i.e. 

N N+1 

(47) £C = i; £x = i. 



n z =0 u=l 



3.4 The permittivity 



While determining dynamic permittivity of crystalline film, the formula by 
Dzyaloshinski and Pitaevski [20] may also be used in the same form (26) in 
which was used for calculating permittivity of corresponding bulk structures, 
with difference that in this case permittivity depends on a film layer n z , i.e: 



(48) 



■(«) = ! 



irS\X\ 
ih 



[G nz {u) + GnA-^)] 



Since starting Hamiltonian was taken in harmonic approximation, ignoring 
(small) members of exciton-phonon interaction [1,6,15], permittivity tensor 
has the real elements only. It is led to that all elements of permittivity tensor 
in one crystalline plane parallel to boundary planes are mutually equivalent, 
i.e. they depend on plane position (n z £ [0,iV]). 

Substituting expression for Green's functions (43) to (46), we obtain ex- 
pression for the relative dynamic permittivity tensor elements in direction 
normal to boundary planes in form: 



(49) 



S 



N+1 



. 1 = 1 -f E E 



v=l s=+, 



where: g± = t|/| — \p\ ~ J~ xy , and after arranging the expression we finally 
follow: 

Qu ~ \p\ + Fxy 



(50) 



£n z (w) 



AT+1 



\f\ 2 -(Qu-\ P \+T xy y 



Figures 5 a-c show a dependence of dynamical permittivity (e) on reduced 
relative energy, i.e. the frequency of external electromagnetic field (/ = 
huj/\X\) for four-layered monomolecular film. Dependence was calculated for 
plane center of the Brilouin's zone (k x = k y = 0), but individually per atom 
planes (parallel boundary areas) of crystalline film, that is for n z = 0, and 
n z = 1,2,3,4. Each graph shows number and position of resonating peaks. 
The resonating peaks on frequency dependence of dynamic permittivity are 
the positions — resonating frequencies where permittivity diverges to ±oo. 
These are also the energies (wavelengths) of electromagnetic radiation which 
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model crystal in given place practically "swallows", i.e. these energies are 
absolutely absorbed there. 

(a) (b) (c) 




«z=0,4 B2=l,3 nfl 



Figure 5. Dynamic permittivity of ideal films 

A final number of resonating peaks is present, because k z component of 
exciton wave vector is discrete. Only for certain values of k z component for 
wave vector (with given values k x and k y ), resonating and radiation absorp- 
tion may occur. Therefore, at least three and maximum five peaks occur - 
number is equal to number of permitted states along axis of translation sym- 
metry breakage, in this case along z-axis. Due to the symmetry of model, 
distribution of peaks in complementary planes is the same: n z = and 
n z = 4 (two boundary planes), as well as for n z = 1 and n z = 3 (first two 
inner planes). Only n z = 2 is different, since for the odd total number of 
atomic planes the middle plane has no complementary plane! 

This result may be explained by experimental facts regarding resonating 
optical peaks in similar molecular layered nanostructures. In papers [23—25] 
this was evidenced in perylene chemical compounds (PTCDA, PTCS, PTFE 
and PBI) and explained by resonating effects at specific unoccupied levels. 
These effects are manifested by narrow optic absorption in close infra red 
band. In comparison with our results, which are concerning deeper infra red 
band in electromagnetic radiation, it may be concluded that these differences 
are effect of differences in crystalline (chemical and physical) structure as 
well as in the model (real and ideal) of samples investigated. 
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4 Conclusion 

This paper describes original application of Paulian Green's function method 
onto theoretical studies of optical properties of molecular crystals. Using 
Boson representation, microscopic (dispersion law and exciton state distri- 
bution) and macroscopic (dynamic permittivity) properties of these crystals 
may be successfully described. 

Strictly following defined procedure, this paper shows that this method 
may be adapted and successfully applied to the study of dielectric properties 
in structures with disturbed spatial-translational symmetry, such as ultra 
thin films. Functioning of newly developed approach is illustrated through 
the problem of finding energy spectra and exciton states, as well as for 
determining relative permittivity of ideal monomolecular film. 

In presence of two parallel boundaries in the system, energy spectrum is 
determined and possible exciton states were found. Important differences, 
comparing with unlimited crystalline structures, have been observed. En- 
ergy spectrum of excitons in monomolecular films is explicitly discrete, and 
the number of discrete levels is equal to the number of atomic planes (in- 
cluding boundary areas) along the axis of the spatial limitations of ultra thin 
film. In bulk sample, a single zone where excitons "take" all possible energy 
values exists. All discrete levels are located within bulk boundaries, and the 
difference in zone width depends strictly (and inversely) on film thickness. 

Comparing with bulk structures, where excitons may be found at any 
place with equal probability, in monomolecular film structures probability 
of finding exciton strongly depends on film thickness. 

In exciton systems of monomolecular crystal bulk, where relative dynamic 
permittivity depends on frequency, continuous absorption zone exists in cer- 
tain range of external radiation energy. In monomolecular film-structures 
resonating peaks exist with precisely determined energies, i.e. resonating 
frequencies. Number of these peaks depends on position of atomic plane 
(regarding boundary planes of the film) for which permittivity is being cal- 
culated: it decreases with the depth of ultra thin film. 

Differences between properties of observed film and corresponding bulk 
structures drastically decrease as the thickness of film is higher. All this 
implies to the existence, and is a consequence, of the quantum size effects. 

Method of Green's differential functions, adapted on described manner 
should be applied further to the study of behavior and properties of more 
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realistic quantum structures, for instance ultra thin films with perturbed 
boundary conditions. 
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